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Abstract 

Substrate and nanotube polarization are shown to change quah- 
tatively a nanotube bandstructure. The effect is studied in a hnear 
approximation in an external potential which causes the changes. A 
work function difference between the nanotube and gold surface is es- 
timated to be large enough to break the band symmetry and lift a 
degeneracy of a lowest but one subband of a metallic nanotube. This 
subband splitting for [10,10] nanotube is about 50 meV in absence of 
other external potential. 

1 Introduction 

Since discovery of carbon nanotubes in 1991 P, a deep physics of these one- 
dimensional nanoscale objects has been demonstrated. Fundamental prop- 
erties of the nanotubes have been studied in view of possible applications in 
electronics and other devices A detailed theoretical description for elec- 
tronic structure of ideal single-wall nanotubes (SWNTs) was obtained as 
well as for an effect of various defects and disorder on the SWNT electronic 
properties ( e.g., Ref.0). However, the nanotube systems under experimen- 
tal study often deviate from a model picture. In this letter we investigate 
one aspect of a real system: the modification of the electronic properties 
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of the SWNT deposited on a substrate. One expects that symmetry of the 
nanotube at the substrate will be lower than the symmetry of the nanotube 
itself in vacuum. 

A description of a breaking of the symmetry of SWNT bandstructure 
due to a charge transfer (or charge injection) between the nanotube and 
the substrate (or contacts) and calculation of a polarization of the substrate 
and the nanotube, which follows to the charge transfer, are the goals of our 
study. Effects of splitting, mixing, and/or anti-crossing of the nanotube sub- 
bands, that are caused by the depolarization of the electron charge density, 
have been almost neglected in literature before. We use term "depolariza- 
tion" for a number of phenomena including a transverse shift of the electron 
charge density from its equilibrium distribution profile (effects due to an ax- 
ial/longitudinal depolarization were discussed elsewhere|ll El E])- We will 
show that the transverse depolarization results in qualitative changes of the 
nanotube density of states (DOS) near van Hove singularities. In particu- 
lar, we predict the splitting of a doublet state ;7J to be likely observable as a 
function of the injected/induced charge density of the SWNT. We will dis- 
cuss that in Sec 12. II In Sec J2.2l we will calculate this injected/induced charge 
density in a self consistent way. 

The depolarization and intrasubband splitting will be studied for a typical 
experimental situation: a single SWNT lies on a conductive substrate or 
separated from the conductor by a thin insulating layer representing an oxide 
on the surface of a metal. We assume that the nanotube is connected to 
electron reservoirs, which may be the leads or the conductor substrate itself. 
A transverse external electric field and/or a work function difference between 
the SWNT and the substrate/contact induce non-zero electron/hole charge 
density in the nanotube. This extra charge density polarizes the substrate, 
which breaks the axial symmetry of the nanotube. This effect is much larger 
than an electronic structure perturbation caused by the lattice distortion 
which may happen due to a van der Waals attraction to the substrate jHj. We 
will demonstrate that a direct action of the uniform external electric field is 
of minor importance as compared to the nonuniform field of surface charges 
on the substrate. We will discuss a modification of our theory of the subband 
splitting for a case of purely insulating substrate in the last section. 
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2 Perturbation Theory for Bandstructure Mod- 
ification 



2.1 Splitting of SWNT subband due to transverse de- 
polarization 

To calculate the splitting and shift of the electron energy levels one needs to 
know matrix elements of the perturbation potential between corresponding 
wave-functions. In our case, the perturbation is a selfconsistent Coulomb 
potential (operator in Heisenberg representation) which describes the inter- 
action between the probe electron and the extra charge density on the SWNT 
and the polarization charge density on the substrate surface: 



Both the probe electron and the nanotube surface charge are taken on a 
cylinder of a radius R. Then z and a are the electron coordinates in the 
cylindrical coordinate system, a is the surface charge density. It does not 
depend on the coordinate Z along the nanotube because we assume the 
translational invariance of the problem for clarity of derivation. Although, 
the theory can be easily extended for the case of slow variation of a along the 
axis. We will show later that one can drop dependence of a on the angle (3 
along the circumferential direction in approximation of a linear response (in 
higher orders of perturbation theory a direct transverse polarization must be 
taken into account 0). cr* is an image charge density which is equal to —a 
for the metallic substrate. 

The first term of Eq.(0) is the interaction with the charge density on the 
nanotube, which coincides with the Hartree term for the SWNT in vacuum 
(without charge injection). The second term in Eq.(^ has also a simple 
physical meaning: this is the energy of interaction of the electron with the 
image charge. The separation between the SWNT axis and the surface of 
the conductor is h. In case of the metallic substrate it is about the nanotube 
radius, R, plus the van der Waals distance for graphite: h ^ R + 0.34 nm. 



V = e 
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{z~ ZY + {Rcosa - Rcospy + {Rsina - 2h ~ Rsin py 



The matrix element of the Coulomb operator is calculated with the 
wave-functions of a tight-binding (TB) Hamiltonian. We use envelope wave- 
functions, obtained similarly to Ref.^^. This approach has been widely used 
in the literature, so we skip details and give the wave-functions in the one- 
band scheme (tt electrons only) in the form: 

\^ljm,k,d = + Cc„fc|5))e^'^e^"^", (2) 

here index m labels subbands of the SWNT electronic structure, k is a lon- 
gitudinal momentum, these two are good quantum numbers (discrete and 
continuum, respectively) for an ideal, long enough nanotube; ( = ±1 is a 
pseudospin. (A pseudospinor vector is formed by a two-component wave- 
function amplitude defined for two atoms in a graphite unit cell, A and B). 
Coordinate along the tube is z, and a is the angle along the nanotube cir- 
cumference. 

We assume that our potential is smooth at the scale of the single unit 
cell (0.25 nm). Then, one may neglect transitions with the pseudospin flip 
(transitions between sublattices). With use of the orthogonality relation 
between the spinor components, it yields: 

{m\V\m) = AireRalog ^"^^ , m = n (4) 

where a, the surface charge density, has to be defined later in a self consistent 
way. 

The Eqs. fl3l4j) are obtained by a direct Fourier transformation of (Q) and 
describe the energy level shift when m = n and the mixing of different 
subbands aX m ^ n. The most interesting term with n = —m is the mixing 
between the degenerate electron states within the same subband. By solving 
a secular equation for the intrasubband mixing of the electron doublet we 
obtain the splitting of the van Hove singularity at the subband edge (Fig|21). 
The new subband energy separation reads as: 

Let us now calculate the injected/induced charge density a which will allow 
us a numerical estimation for the 5Em splitting. 
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2.2 Charge injection due to the Fermi level shift 

The Eqs.(jSHSl) written for the given charge density a which has to be 
derived in this section. When the SWNT is placed in a real device, one must 
consider the work function difference between the nanotube and the contact 
or the conducting substrate and/or the external potential which may be 
applied to the SWNT. The potential shifts the Fermi level in the SWNTj^]. 
As a result of this the positive/negative charge is injected into the nanotube: 



a 



2txR 



iy{E)dE, (6) 







here z^(-E') is a bare one-dimensional DOS (independent of a in a linear re- 
sponse theory); n = AW — e<^^* — e(p^"''^{a) is the shift of the electrochemical 
potential of the SWNT (with respect to a charge neutrahty level E = 0) 
which depends on the work function difference, AW, on the external poten- 
tial, (y?^*, applied between the nanotube and the reservoir and on the potential 
^md induced by the charge density of the nanotube, a. This last term is pro- 
portional to the intrasubband term (m = n) of the Coulomb interaction given 
by the Eq.®. 

This selfconsistent equation for a, is readily solved analytically if the 
electrochemical potential is below the second subband edge. We follow Ref.jl] 
in derivation of a: the induced potential is obtained by direct integration of 
the charge density along the SWNT as in IH IH] . As it is shown in Figure 
for a metallic SWNT, the charge is a product of the constant DOS, Cq, and 
the electrochemical potential, fi. Then the solution of Eq.® is as follows: 

AW-eip''^ 

= 7 TV, (7) 

271 Re [2 log{2h / R) +Cq^) 

and for a semiconductor SWNT, which has the DOS oc CQEe{E-A) /^E'^ - A^, 
the charge is as follows: 



y (aW^)' (41og^(2/^/i?) - Cq-") + Cq' - 2 \og{2h/R) 

(8) 
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Here 9{x) is the Heaviside unit step function and A is 1/2 of the energy gap. 
We introduced a quantum capacitance of the SWNT following Ref.jl]: 



Se^ 

which is the one-dimensional analog of the quantum capacitance proposed 
for a two-dimensional electron gas system by Luryi^I]. Here b ~ 1.4A is the 
interatomic distance, 7 ~ 2.7 eV is the hopping integral for the graphite-like 
systems. We notice that despite the az, as given by Eq.(jH)), comes from a 
massive subband (in contrast to a a, as in Eq.Q where the lowest subband is 
massless, see in FigHj), the linear dependence of <7a/z on ip^^ preserves as long 
as the potential yj^* is large enough. This reflects the fact that a classical 
one-dimensional charge density is a linear function of a classical electrostatic 
potential 



3 Results and Discussion 

In the last section we obtained the selfconsistent expression for the surface 
charge density as a function of the external potential and the work function 
difference which may be considered as a built-in potential. Substituting 
Eq.((Zj) into Eq.(jn)) we obtain the splitting of the degenerate subbands | ±m) 
of the metallic SWNT (when the Fermi level is within the first subband) as 
follows: 

The splitting decreases with m exponentially, hence, the effect is likely ob- 
servable for the lowest degenerate subband. Then, for the following param- 
eters: the SWNT radius R ~ 6.7 A, the distance to the metal substrate 
h = 10.1 A, and the quantum capacitance Cq^ ~ O.GQjlj, we obtain a numer- 
ical estimate for the subband splitting SEi ~ 0.15 {AW — ey?^*). Experimen- 
tal data for the work function of SWNTs scatters from 4.9 to 5.05 eV jT^E]- 
For the SWNT on the gold substrate we use as an estimate AW ~ 0.3 eV. 
In absence of the external potential, this work function difference results in 
~ 46 meV gap between two spht peaks of the density of states (FiglSJ, which 
is larger that kT at room temperature. We also calculated the contribution 
of all other subbands, which is negligible in the splitting, but it shifts the 
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doublet as a whole. As a result, two new peaks in the Fig|21 appear not 
symmetrical with respect to the original DOS singularity. 

The splitting of ±m doublet is an analogue of a degenerate level Stark 
effect for the nanotube in a multipole potential of the image charge. The 
lower subband has x symmetry and the upper subband has y symmetry 
(with corresponding wave-functions |x) = l/^/2{\ + m) + \ — m) and \y) = 
l/^/2{\ + m) — I — m)) because of an attraction energy of the electron to its 
image charge is lower for the second combination. 

We predict a similar effect for the semiconductor nanotube, although, the 
total external potential causing the charge density injection must be larger 
than one half of the gap in this case. As we study in this paper only the 
effect which is linear in the external potential, all high order terms in Eq.® 
have to be discarded. 

3.1 Dipole polarization correction 

The charge injection in the nanotube may be readily achieved by applying 
external electric field. One may naively argue that the external field itself 
can break the bandstructure symmetry and result in some level splitting. 
Although, it is correct statement in general, the direct splitting of the SWNT 
orbital doublet ±m by the uniform electric field is forbidden by symmetry. 
These degenerate states do not mix together due to the selection rules of the 
problem. The matrix element for an intrasubband splitting in the uniform 
external field S^t equals zero by parity: {m\e£xty \ — m) = 0. 

In order to calculate the subband splitting in this case we have to compute 
the charge injection, which is proportional to the applied field. The potential 
is equal to y?^* = Sxth, where h is the distance between the axis of the tube 
and the metal surface, which has to be substituted into Eq. ljlUj) . 

In Sec J2.1l we assumed that the charge density a has no dependence on the 
angular coordinate f3 along the nanotube circumference. This is an accurate 
approximation since a dipole (and higher multipole) component of the a is 
small as compared to what is given by Eqs. (|7l8| ). Let us prove this assumption 
for the dipole polarization of the SWNT. 

The non-uniform external potential causes a deviation of the surface den- 
sity from the uniform equilibrium value, a, which is given by following ex- 
pression: 

5aiP) = eY. ^^'-J^^^f^'\ mm, (11) 
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where fi are the occupation numbers, the matrix element (^|V^|j) is given 
by Eqs. (j3l4p . Ei are the energies of subbands and are corresponding 
wave-functions. 

We define the nonuniform dipole part of the charge density of a SWNT 
as 6ai = J sin Pa{P)dp. Then, the dipole component of the surface charge 



is as follows: 

Let us remind that according to Eq.Q {i\V\i ± 1) = —i8TTR'^ea/{2h). 

In case of the electrochemical potential equals zero (no charge in the nan- 
otube), the transverse polarization includes transitions from the valence to 
conduction band only (t>|V^|c) (the details of the calculation are presented 
elsewhere j9j ) . Here, we study an extra component of the polarization, which 
is due to the induced charge density. Thus, we need to consider only tran- 
sitions from the levels above the charge neutrality level, E = 0, and below 
the Fermi level, E = Ep, (which is the shaded area in Fig^. Hence, the 
dipole polarization is proportional to the net charge density a, and the dipole 
charge density of the armchair SWNT is given by the following expression: 



(2nRaAf R , 2h 
—- log — 

327r e h ^ R 



5a,= ^-Q^—^ l^iog^ oc^^, (13) 



where Cq ~ 3.2 is the dimensionless quantum capacitance. We single out 
term 27rRaA, which is the specific one-dimensional charge density of the 
SWNT and proportional to the external potential and, thus, to the external 
field. 

The Eq. ()13j) shows that the effect of the transverse polarization on the 
bandstructure is quadratic in the external field, in good agreement with a 
plain dielectric response theory [0], while the effect of the image charge is 
linear in Sxt- Thus, the degenerate level splitting due to the dipole component 
of the polarization will be less important than the splitting due to a uniform 
component: ctq = lo^ cr{P)dP, at least, in a weak field regime discussed in 
the article. This proves post factum our assumption of a to be independent 
oi/3. 
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3.2 Depolarization at the insulator substrate 

For the sake of completeness we present here also a modification of our theory 
to the case of a dielectric substrate. In this case the screening of the charge 
density in the nanotube is weaker. It results from (i) underscreening of 
the Coulomb interaction between the nanotube carriers and (ii) lower charge 
density induced in the substrate. The second factor can be taken into account 
by substituting an effective image charge density a* = in the second 

term of Eq.(IT)), where e is the dielectric function of the substrate (in case of 
highly conductive substrate it equals — oo), instead of the bare image charge 
density —a. This results in substituting a* in Eqs.(j2HHl) where appropriate. 

Now, the fields of the image charge and the charge in the SWNT do 
not cancel each other in contrast to the case of the metallic substrate. As a 
result, the underscreening of the Coulomb interaction happens. This modifies 
the equations for the energy level shift ( intmsuhhand matrix elements as in 
Eq-0) and, thus, the electrochemical potential shift. One must substitute 
\og{2h/R) term everywhere by \og{2h/R) + 2/(e — l)log(L/i?) where L is 
the length of the nanotube (or distance between metal leads to it). This 
expression diverges with the length of the nanotube which reflects the one- 
dimensional character of the Coulomb interaction. These changes have to be 
made through Eqs.fffHTUj). 

The first term of Eq.(PQ) does not appear in the calculation of the inter- 
subband matrix elements as in Eqs. (|3l5|) . Hence, no additional correction is 
required in the equations (fTTHT^ of the last section. 

We assumed in this paper that the perturbation theory in a linear ap- 
proximation in fi (or equivalently in a) is applicable. Restrictions which may 
follow from this assumption are as follows. The external potential has to be 
small. We neglect here the dipole term in the induced charge density (and 
higher multipoles as well). It is equivalent to a weak intersubband mixing 
which assumption may not hold for wide nanotubes or strong external fields. 
The effect of the strong field on the bandstructure is discussed elsewhere 0. 
In this paper we used the Eq.ijHj) for the equilibrium charge density in the 
SWNT. One may consider transport devices on equal basis, as long as the 
charge of the nanotube is still given by the quasi-equilibrium charge density. 
However, for non-zero current flowing through the nanotube, one must use 
an expression for the charge density which differs from Eq. (jH)) (to be discussed 
elsewhere |12j). 
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4 Conclusions 



In summary, we have developed a microscopic quantum mechanical theory 
for a charge transfer between a SWNT and a conductive substrate (and/or 
metallic leads). This charge injection results from a natural work function 
difference between the nanotube and the substrate or/and from an external 
potential applied between those. A surface charge density of the SWNT is 
calculated selfconsistently within an envelope function formalism of tight- 
binding approximation. 

We demonstrated for the first time that the influence of this charge trans- 
fer on the electronic structure of the SWNT is not negligible for typical 
material parameters of the problem. Because of the brcak;ing of the axial 
symmetry of the system, the SWNT DOS changes qualitatively: degenerate 
subbands ±m, where m 7^ 0,n, split. It has a simple physical interpretation 
— the electrons with x and y polarizations are no longer equivalent as their 
attraction to the substrate is different. This effect can be related to a de- 
generate level Stark effect with an appropriate choice of external field of the 
image charge. The gap between the new x and y subbands is constant in 
/c-space (for the external field which is uniform along the tube) so it shows 
up dramatically at the subband edge. The van Hove singularity splits, and 
the distance between two peaks of the DOS is about 46 meV for the [10,10] 
armchair SWNT on the gold sul:)stratc. 

We obtained analytical expressions for matrix elements of the image 
charge field, which yields the mixing of different subbands and can be used 
to describe the level anticrossing (to be discussed elsewhere). Same ma- 
trix elements enter the expression for the multipolc polarizabilites of the 
SWNT. We estimated a major contribution to the dipole polarizability of 
the armchair SWNT, which comes from intraband transitions for non-zero 
charge injection. The analytical expression for the dipole component of the 
surface charge density is shown to be proportional to the square of the ex- 
ternal potential and, hence, appears in the second order of the perturbation 
theory which corroborates post factum our assumption of uniformity of the 
induced/injected charge along the SWNT equator. 

We show that the modification of our theory to the case of semiconductor 
substrate is straightforward. The analytical expressions for the van Hove 
singularity splitting and induced charge density are obtained. 
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Figure 1: Schematic density of states (DOS) of a metallic SWNT. First 
(massless) subband contributes to a constant DOS at the E = 0. When the 
Fermi level, Ep, is lower than the second (massive) subband edge (which 
corresponds to the first peak of DOS), an injected/induced charge is propor- 
tional to the shaded area and is a linear function of Ep. 
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Figure 2: Density of states (DOS) of [10,10] armchair nanotube in vicinity 
of first van Hove singularity (black color). Charge injection in the NT due 
to work function difference (see the text) results in a splitting of a doublet, 
which is clearly seen as compared to bare DOS of neutral NT (light gray 
color) . 
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Figure 3: Image for table of contents: Zigzag [17,0] nanotube on a surface of 
Si02 substrate. Geometry has been relaxed with Molecular Mechanics. 
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